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LONG  TERM  GOALS 

The  long  term  goal  of  this  project  is  to  develop  efficient  inversion  algorithms  for  successful 
geoacoustic  parameter  estimation,  inversion  for  sound-speed  in  the  water-column,  and  source 
localization,  exploiting  (fully  or  partially)  the  physics  of  the  propagation  medium.  Algorithms  are 
designed  for  inversion  via  the  extraction  features  of  the  acoustic  field. 

OBJECTIVES 

•  Achieve  accurate  and  computationally  efficient  inversion  for  propagation  medium  parameters  and 
source  localization  by  designing  estimation  schemes  that  combine  acoustic  field  and  statistical 
modeling. 

•  Develop  methods  for  passive  localization  and  inversion  of  environmental  parameters  that  select 
features  of  propagation  that  are  essential  to  model  for  accurate  inversion. 

•  Implement  Bayesian  filtering  methods  that  provide  dynamic  and  efficient  solutions  for  the  first  two 
objectives. 

APPROACH 

Continuing  efforts  from  previous  years,  we  worked  with  Bayesian  approaches  applied  to  sound  signals 
for  the  extraction  of  acoustic  features  using  a  combination  of  physics  and  statistical  signal  processing. 

One  of  the  topics  approached  this  past  year  was  source  localization  and  water  column  sound  speed 
estimation  using  arrival  time  estimation  for  propagation  in  multipath  environments  with  sequential 
Monte  Carlo  methods,  namely,  particle  filtering.  Some  of  our  results  have  been  presented  in  [1,  2,  3] 
and  [4];  the  latter  paper  is  in  press.  The  method  is  an  extension  of  previous  approaches  developed  by 
the  PI,  now  using  more  advanced  techniques.  The  initial  goal  is  to  estimate  accurately  arrival  times 
and  corresponding  amplitudes  of  sound  paths  in  shallow  water  environments.  Arrival  times  as  well  as 
amplitudes  of  corresponding  paths  provide  critical  information  on  the  geometry  and  environment  in 
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which  sound  propagates.  Even  small  discrepancies  in  arrival  time  and  amplitude  estimates  can  lead  to 
significant  errors  in  inferences  on  inversion.  Accurate  estimation,  however,  linked  with  physical 
models  for  sound  propagation  in  the  ocean,  can  lead  to  accurate  source  localization  and  inversion  for 
properties  of  the  propagation  medium.  As  will  be  seen  below,  in  addition  to  simple  point  estimates, 
our  method  facilitates  the  computation  of  posterior  probability  densities  (PPDs)  for  arrival  times, 
which  can  be  used  for  uncertainty  quantification  in  inversion  (source  localization  and  depth  and  sound 
speed  estimation  in  the  water  column;  sound  speed  estimation  in  the  seabed  has  also  been  performed). 
These  distributions  can  be  multimodal,  a  fact  that  is  not  often  revealed  when  different  methods  are 
used.  It  has  also  been  shown  in  [1]  that  amplitude  PPDs  can  be  obtained  in  a  straightforward  manner; 
those  can  be  used  for  seabed  attenuation  estimation. 

More  work  was  conducted  on  identifying  arrival  times  and  amplitudes  of  distinct  frequencies  in 
dispersion  curve  analysis  (within  a  single  mode  or  across  different  modes).  Those  characteristics 
provide  significant  information  on  properties  of  the  propagation  medium  and  source  location.  The 
importance  of  modal  arrival  times  and  amplitudes  in  geoacoustic  inversion  and  source  localization 
using  dispersion  curves  has  been  extensively  discussed  in  [5,  6,  7,  8], 

WORK  COMPLETED 

Previously  developed  approaches  were  improved  to  provide  more  accurate  estimates  of  arrival  times 
and,  subsequently,  of  source  location,  bathymetry,  and  the  water-column  sound  speed  profile. 
Specifically,  in  addition  to  conventional,  forward  particle  filters,  smoothing,  backward  mechanisms 
were  employed,  that  reduced  uncertainty.  The  arrival  time  tracking  method  was  applied  to  synthetic 
and  real  data  for  source  localization,  and  bathymetry  and  sound  speed  estimation  in  the  water  column 
and  the  seabed  layer.  Results  were  obtained  by  using  the  particle  filtering  results  as  input  to  the 
inverse  problem;  more  results  will  be  obtained  this  year  by  further  improving  and  testing  our  approach. 
Passive  fathometer  processing  with  particle  filters  was  also  improved,  with  the  developed  methodology 
applied  to  seismic  tremor  tracking  as  well  [9],  Additionally,  the  particle  filtering  method  was  applied 
to  synthetic  long-range  data  for  dispersion  curve  extraction  with  results  significantly  improved  in  terms 
of  resolution  from  those  presented  in  [10]. 

RESULTS 

Results  that  we  have  previously  reported  illustrated  how  we  can  extract  arrival  time  estimates  from 
sound  signals  in  the  ocean  that  are  more  accurate  than  those  computed  with  conventional  methods. 
However,  perfonning  inversion  with  arrival  time  estimates  obtained  with  our  method,  we  noticed  that, 
in  some  cases,  uncertainty  in  estimates  can  still  be  significant. 

To  improve  our  results,  we  further  extended  our  arrival  time  tracking  method  to  obtain  tighter  PPDs  on 
arrival  times.  Reducing  the  variance  in  the  arrival  time  PPDs  is  directly  related  to  limiting  the  variance 
in  source  location  and  environmental  parameter  estimates.  We  observed  that  a  portion  of  the 
uncertainty  in  our  arrival  time  estimation  is  a  result  of  limited  prior  infonnation  in  the  initial  states  of 
the  filtering  process.  To  mitigate  this  problem,  we  employed  backward-moving  filters  that  improved 
results  on  initial  states  by  moving  backwards  after  estimation  at  all  states/receivers  was  completed 
following  the  approach  of  [1 1].  We  extended  our  work  in  this  direction  during  the  last  year.  Figure  1 
shows  results  from  an  extensive  perfonnance  evaluation  of  the  smoothing  process  as  compared  to  a 
standard  forward  only-filter. 
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Figure  1:  Demonstration  of  the  advantage  of  a  PF  with  smoothing  over  a  simple  PF  for  three  paths 
at  an  SNR  of  10  dB:  (a)  PF-MAP  RMS  errors  for  arrivals  times  from  a  forward  PF  for  three  paths 
(circles  for  D,  stars  for  SR,  and  triangles  for  BR)  and  a  forward-backward  filter  for  the  same  paths 
(asterisks,  diamonds,  squares);  (b)  RMS  errors  for  arrivals  times  for  a  forward-backward  filter  for 

the  three  paths  (*)  and  GS-MAP  results. 


Specifically,  the  figure  shows  Root  Mean  Squared  (RMS)  errors  for  Maximum  a  Posteriori  (MAP) 
arrival  time  estimates  for  three  paths  for  (a)  a  forward  filter  and  a  forward-backward  filter,  and  for  (b) 
the  same  forward-backward  filter  and  the  Gibbs  Sampling  (GS)-MAP  processor  presented  in  [12]. 
Results  show  that  the  forward-backward  filter  is  superior  to  both  the  forward-only  filter  and  the 
conventional  processor. 

To  further  illustrate  the  power  of  the  process,  we  show  in  Figures  2  and  3  source  localization  and 
bathymetry  estimation  PPD  results  from  applying  the  method  of  [12]  and  the  new  method  to  the  Haro 
Strait  data.  Figure  2  shows  PPDs  obtained  with  the  GS  method  of  [12].  Figure  3  presents  results 
obtained  with  the  new  filter.  The  reduction  of  the  uncertainty  in  the  PPDs,  evident  by  comparing 
Figures  2  and  3,  is  significant. 

Additional  work  was  carried  out  by  applying  particle  filters  to  real  data  for  sound  speed  estimation. 
Specifically,  we  worked  with  data  from  the  SW06  experiment,  collected  at  the  16-element  MPL-VLA1 
array  (because  of  low  SNR  at  the  15th  and  16th  phones,  data  at  only  the  14  lower  phones  were  used) 
[2],  The  source  signal  was  a  linear  frequency  modulated  pulse  with  frequencies  between  100  and  900 
Hz.  Figure  4  illustrates  estimation  of  the  sound  speed  profile  for  the  water  column  using  the  ray  path 
arrival  time  PPDs  such  as  those  shown  in  [2],  Specifically,  Figure  4  shows  the  PPD  of  the  coefficient 
of  the  first  EOF,  which  is  in  agreement  with  results  presented  in  [13].  Note  that  the  PPD  is  bimodal, 
which  is  not  evident  in  results  previously  presented  in  work  with  the  same  data.  The  resulting  sound 
speed  profile  was  similar  to  the  profile  shown  in  [14]. 
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Figure  2:  PPDs  of  (a)  source  range  and  (b)  depth  and  (c)  water  column  depth  estimates 

with  the  Gibbs  Sampling  method  of  [13]. 


The  Pi’s  particle  filtering  method  for  dispersion  tracking  was  improved  working  along  similar  lines. 
The  new  method  was  applied  to  synthetic  long-range  data  for  dispersion  curve  extraction  with  results 
significantly  improved  from  those  presented  in  [10].  These  results  are  shown  in  Figure  4,  where  MAP 
estimates  are  superimposed  on  a  spectrogram.  The  method  detects  modal  trajectories  that  exit  at  low 
frequencies  and  new  that  appear  at  high  frequencies.  The  continuous  or  interrupted  nature  of  the  curve 
captures  the  degree  of  reliability  of  specific  modal  behaviors  in  inversion. 
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Figure  3:  PPDs  of  (a)  source  range  and  (b)  depth  and  (c)  water  column  depth  estimates 

with  particle  filtering. 
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Figure  4:  Sound  speed  estimation  for  the  SW06  experiment  using  arrival  times  extracted  by  the 
particle  filter  and  used  as  input  to  the  inverse  problem:  the  figure  shows  the  PPD  of  the  coefficient 
of  the  first  EOF.  The  results  are  in  agreement  with  those  presented  in  [13,14]. 


IMPACT 

The  significance  of  accurate  arrival  time  estimation  in  source  localization  and  water-column  sound 
speed  estimation  has  been  studied  with  the  goal  of  producing  accurate  parameter  estimates  and 
measures  of  the  uncertainty  in  the  estimation  process.  The  reliability  of  arrival  time  estimation 
methods  is  intimately  tied  to  the  quality  of  the  inversion  process.  The  new  approaches  facilitate  arrival 
time  extraction  and  the  association  between  paths/modes  and  detected  arrivals  and  also  produce  PPDs 
of  arrival  times,  modal  frequencies,  and  corresponding  amplitudes. 

RELATED  PROJECTS 

The  PI  is  collaborating  with  Drs.  Yardim  and  Gerstoft  on  sequential  filtering  in  ocean  acoustics.  The 
PI  is  also  collaborating  with  Dr.  Leon  Cohen  on  comparing  numerical  and  analytical  descriptions  of 
dispersion  in  ocean  environments  and  is  involved  in  discussions  with  Dr.  Ross  Chapman  on  inversion 
for  attenuation. 
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Figure  5:  Dispersion  curves  estimated  with  particle  filtering  from  a  synthetic  environment 
simulating  the  environment  of  the  Gulf  of  Mexico  experiment. 
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